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We study and solved Schwinger-Dayson equation for massive quarks in the quark 
model. Interaction potential is choose as sum oscillator and Coulomb terms. The cal- 
culation show us that negative energy of quark exist under m <C (4V D 3 /3) 1 ^ 3 i.e. when 
mass quark is less that parameters of scale. Us are offered to use asymptotical behaviour 
at the solutions Schwinger-Dayson equation selection physical and non-physical condi- 
tions. The quark condensate under different values mass quarks are calculated. The scale 
parameters (4V 3 /3) 1/3 = 520 MeV. 
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I. INTRODUCTION. 



For description of decays of particles to us needed will know a behaviour of wave functions and the 
spectra of masses of mesons and baryons. In literature there are several general investigation ways: quark 
model, lattice QCD, bag model, QCD sum rules, quenched QCD and other phenomenological models and 
approaches. They have their own advantages and defect. 

The most real approach is quark model based on Schwinger-Dyson (SD) and Bethe-Salpeter (BS) 
equations. This quark model proposed by the Orsay group [El, perfection by Pervushin Q and Bicudo- 
Ribeiro [Hj. The main purpose the works [|l]-[| was offer a method define wave functions and obtain 
spectrum masses of mesons and hadrons. Next step development of quark model is calculate meson 
decays [Q] and improve upon |^| . Recently in the paper || presented the development of the quark model 
and they have point that model a pretender on the description hadron physics. 

Thereby we have quark potential model which based on SD and BS equations. The definition of the 
form of quark-antiquark potential is one of the actual problem of the meson spectroscopy. 

In the present paper we investigate quark mass function based on SD equation. We consider quark in 
hadron with interaction potential chosen as sum Coulomb and oscillator terms. We has solve SD mass 
gap equation and the energy of quarks are calculated in this approach. Using asymptotic behavior when 
p — * oo we can select solution on physical and non-physical stages of quarks. By means of solutions SD 
equation we have calculate < qq > quark condensate under different masses of quarks. 

Let us consider the effective covariant relativistic action for singlet channel 



S eff= I d 4 xq(x)(ifl - m )q(x) - - / d' l xd' l yq l 3 2 {y)q ai (x) 



vi \ X + V\ 
K\{x-y\—^-) 



9/3! {x)q a2 (y) 



(1) 



with interaction (kernel) 



Ki(x - y\X) au/3l/a202 = fi ai01 [V(zj_S(z ■ v)]fl a2/ 3 2 

where rh = (m G i, m onf ) are the bare masses, a or (3 are the short notation for the Dirac and flavour 
index. Quantization axis is chosen 
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where rj^ - single time-axis. From the covariant effective action ([!]) we can obtain SD equation for mass 
operator (E(p^)) depends on transversal momentum p 1 - with nonzero mass current quark (m D ) @j. 



where j) — 77^7^, rj^ is satisfy Markov- Yukava principle (Q), 

p^ = p^ — i]^(r]-p)- momentum, which orthogonal on time axes (p±, r\ = 0), 

and Gs - Green's function of quark 

Gx(q) = l4-X(qT 1 (4) 

Let us note that after integration (||) over " time" component (g M = q ■ rf) we obtained the mass operator 
with perpendicular component only. We can write in Lorentz covariant presentation : 

= # x + E(p ± )S- 2 (p ± ) (5) 

where E - energy of quark, and S - unitary operator which can to present (parametrize) through v function 
as 

S(p±) = exp{-f,v{p±)} (6) 



where p± = p±/\p±\ , \p±\ = \J (p±_) 2 - The one-particle Green function (Q) can be represented in the 
form of the expansion over the states with positive and negative energies 
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G s (g) 



A r ( ' +) (^) 



9 f >-£% ± )+ie <f' + £(g ± )+ie 



(7) 



(8) 



A^_ ) =5( g + )A^ ) (Q)5- 1 (^) 
where Ay_^ (0) are projection operators of positive and negative value, 

A? +) (0) = i(I + tf), 



(9) 



(10) 



A^O) = -(!-#), 



(11) 



here I - one-dimensional 4-matrix. 

Since energy of quark E(Q) ~ m , E(oo) ~ and considering that 5 2 (0) = 1, S' 2 (oo) = —1, so 
propagator of quark in limit p± — > coincide with propagator free quark in rest system 



G(0) 



1+7° 



1-7° 



2 — to — ie p|| + to — ze 
In limit pj_ — > oo propagator G(oo) coincide with propagator massless quark 



G(oo) 
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1 + 72? 



1 - TP 



;>L> cquat' 

is quark. Inserting (U)-(0) into ^3|) and after integration over q v = q° we have get next equation for v 



Let us consider SD equation (M) in rest system 77^ = (1,0,0,0), of bound state, one of the praticles 
quark, 
function 

d 3 p 



1 f drp 

m o sin 2v (p) — p cos 2v(p) = - / V(\p — q\)[cos2v(q)sm2v(p) 

2 J (27r) J 



p ■ q 

sin 2v(q) cos 2v(p)] 

pq 



(12) 



and expression for energy of quark 



E(p) = m cos 2v{p) + p sin 2v{p) + 



d 3 p 



■> / (2tt) 3 Ql)[co52v( 9 ) cos2w(p) - 

p • ^ 

sin2w(g) sin 2v(p)] 

pq 



(13) 



II. SCHWINGER-DYSON MASS GAP EQUATION. 
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In the literature the equations (|12|),(|13[) usually called as mass gap equation. Lets us make v(p) 



tp{p)] variable change in (|l 2|) (jT3|) . So we obtain the following equations 

E(p)sintp(p) = m Q + ^ J -^-^V(\p- q\) sin ip(q) 
1 f aTa _ 



E(p) cos ip(p) = p + 



(14) 



2 J (27T)' 



fV(b- g|)pgcosy>(cj>) 



where p = |p| , p = 
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The SD equation define one-particle energy of flavour quark. The potential phenomenology of spec- 
troscopy of quarkonia Q] in the meaning of the real interaction, uses in general, the sum Coulomb and 
the increasing potential. In the increasing potential is used the lattice calculation for description of the 
heavy quarkonia. The definition of the form qq - potential is one of the actual problem of the meson 
spectroscopy. We chose the potential in the form as sum of oscillator and Coulomb terms. 



V{p) = - : 



47ra s 
p2 



One particle SD equation ( ]14| ) for quark in dimensionless units takes a form 
d 2 ip(p) 2 dip(p) s'm2(p(p) 



p A 



2p sin ip(p) + 2m cos <p(p) 



dp 2 p dp 

+ I dq- {Tli{p,q) cos ip{p) - lZ 2 (p, q) shnp(p)} = 
37T J p 



(15) 



(16) 



where 



E(p) —pcoaip(p) +m smip(p) — —[<f'(p)Y 



COS 2 lfi(p) 



2a s 
' 3tt 



J dq p {^-(p> sin v(p) - ^(p, q) COS (f(p)} 



Ki{p,q) = In 



P +q 
2pq 



In 



siny>(g) 
p + q 



p-q 



m 
r(q) 

- 1 



cos f(q) 



r(q) 



(17) 



r(q) = \/ q 2 + (m ) 2 . The equation ( |l6| ) and expression for self-energy of quarks ( |l7| ) is present in 
dimensionless unit (4V^ ) 3 /3) 1 / 3 = 1. We can define this parameter in physical region after solve BS 
equation, from spectroscopy mesons and baryons. In literature [Q, g], 0] usually use as: (4V r Q 3 /3) 1 / 3 « 
280 -j- 480 MeV. 



III. NUMERICAL METHOD. 



To solve the equation ( fig ) by analytical method is difficult that is way we'll find the solution it by use 
numerical method. Rewrite equation ([l6|) in following type: 



Fi[<p(p)]+F»[<p(p)] = 0, 



(18) 



where 



d 2 ip(p) 2d(p{p) sn\2ip(p) 



dp 2 



p dp 



— 2p sin tp(p) + 2m cos (p(p), 



K[<P(P)] 



2a s 
3tt 



dq- {Ui(p, q) cosip(p) - K 2 (p, q) simp(p)} . 



The problem can be solved step by step scheme: 

i) a s — 0, m = 0, (4V r Q 3 /3) 1 / 3 = 1. The current quark mass equal null end equation ( fl8l ) disregarding 

Coulomb interaction potential is considered. We have chosen interaction potential as increasing 
oscillator potential. 

ii) a s = 0, m ^ 0, (4V^ 3 /3) 1 ^ 3 = 1. The equation ( |l8| ) with provision for current quark mass on the 

intervale ni Q ^ ^-Largemassquark' 



Hi) a s ^ 0, m ^ 0, (4V r 3 /3) 1 / 3 = 1 The equation ( |l8| ) oscillator plus Coulumb term of interaction 
potential and current quark mass not equal zero. 
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Unfortunately to consider the equation ( jig ) for case Hi) without oscillator terms of interaction potential 
(i.e. (4V D 3 /3) 1/3 = 0) to us it was not possible. 



We have derived numerical solutions to the integr-differential equation (18) with boundary conditions 



V(0) = | , ^(oo) = (19) 

using the computational scheme developed in ref . ]7j . This scheme consist the modification of the Newton 
iterations by combining the continuation over a parameter (m Q or a s ) method. This scheme has a 
convenient algorithm for assignment of the initial conditions as functions of the external parameters (m D 
or a s ),due to a special choice of the iteration parameter as a step of the modified Newton iteration. So, 
in this scheme the solutions to problem (p^,(|l9|) are used as an input. 



The numerical solutions to boundary problem (|18|),(19) for several values of m are shown in Figure 
1, under a s — 0.2. The values of free parameters m , (4V^ /3) 1 / 3 and a s belonging to the physical region 
will be defined from the solutions to the BS equation for mesons. 



IV. RESULTS AND CONCLUSION. 

We solve SD equation by using numerical method with choose potential ( |l5|) for massive quarks. In the 
event of massless quarks (m — 0) and without Coulomb term of potential our results coincide with H. 
The SD equation have solution with some nodes. The asymptotical behaviour ip when p — > oo is positive 
value because tp(p — » oo) > for massles quark and tp(p — > oo) = The authors Q study asymptotic 
behaviour when p — > oo for massles quarks and without Coulomb term of interaction potential. 

We consider without nodes solution of SD equation, though needed to note that equation ( |l8| ) has an 
infinite number of solutions. From asymptotic behavior ip we are to select on physical and non-physical 
solutions. Behavior ip under p — > oo following: 

<p(p -> oo) -> ^ (20) 
P 

and it has positive sign, this fact indicates us that the nodes solution of number n = 0, 2, 4, ... (i.e. even 
number) they are physical. The nodes solution of number n = 1,3, 5, ... makes no physical sense because 
asymptotic condition ( p0| ) is violated. 

We have calculate self-energy of quarks under different parameters mass, see Figure 2. Studies have 
reveal us behavior E(p) under m > (4U D 3 /3) 1 / 3 with certain accuracy comply with nonrelativistic ap- 
proximation \J p 1 + m*, authors Q were calculated such results. Presence of negative energy (when 
m = 0) note on Authors || researched a negative energy and adding in the potential of a constant 
energy shift they verified BS equation for mesons and baryons. In fact jlj adding to the potential in 
space coordinates a constant energy shift then the quark energy is shifted by exactly half of that constant 
opposite sign. 

We consider a region when m < (4U Q 3 /3) 1 / 3 . On the fi gure 2 we see that negative energy has an end 
value (under p = 0,m o = 0) E(0) = -6(4V 3 /3) 1/3 , reference to §. With increasing masses an energy 
of quark grows in the negative area and for m > (4U Q 3 /3) 1 / 3 continues to grow in the positive area. On 
the Figure 2 (line (a)) have plot a curved line E(p) = p — in the neighborhood p = 0. 

In our approximation the quark condensate define through the part nonperturbative Green's function 



d 3 p 
(2n) 



< qq >= i2N c tr \ ^[G E (p) - G m (p)} 



J \M) VP + m o 

where Gs(p) = \p — S(p)] -1 and G m (p) = — m)~ 1 are Green's functions of "dressed" and "bare" 
quarks respectively, and tp(p) - the solution of SD equation of quark q. 

In the paper |l[, as we already mentioned earlier, SD equation is considered in the case m = and they 
calculated < qq >= — (140Afeu) 3 in units (4U Q 3 /3) 1 / 3 = 289Meu, we have repeat this result. Moreover 
have calculate behaviour < qq > under different values of masses m and unit parameter (4V r D 3 /3) 1 / 3 , see 
Figure 3. In the event (4V Q 3 /3) 1 / 3 = 520 Mev (line (c) in the Figure 3) our results compare excellent with 
commonly value of quark condensate (-250 Mev) 3 Hj. Completely proper value of scale (4U G 3 /3) 1 / 3 will 
be calculate after solving the equation BS in our approach. 
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FIG. 2. Quark self-energy E(p) is calculated from <p(p) under different quark masses in dimensionless unit. (a)In 
the neigborhood p=0 and a s — behaviour E(p) from expression (jl^) is E(p) =p — -p\ (b) under m — 0,a s = 

repeat results Orsay group; (c)massive quark m = 0.2(4V o 3 /3) 1 '' 3 , without Coulomb term potential a s = 0; (d) 
massive quark m = 5(4V 3 /3) 1//3 , potential as sum oscillator and Coulomb term a a = 0.2; (e) massive quark 
m = 5(41/ 3 /3) 1/3 , without Coulomb term a s = 
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FIG. 3. The quark condensate under different masses quark m = m i +rn 2 (a)under (4V 3 /3) 1/3 = 289 Mev in 
point m = this results complies with Orsay group; (b) under (4V D 3 /3) 1/3 = 368Mev; (c) under (4V; 3 /3) 1/3 = 520 
Mev. 
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